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I. INTRODUCTION 

Neutrino physics becomes one the most interesting 
fields of elementary particle physics, especially in the 
wake of the recent experimental achievements in study- 
ing of solar and reactor neutrinos (see, e.g., Refs. 
One of the most intriguing puzzles in neutrino physics 
in the neutrino oscillations problem. Neutrinos of one 
type can be converted into another type if there is a 
mixing between different neutrino eigenstates. Presently 
three major opportunities for neutrino oscillations are 
discussed in literature. The first one is neutrino flavor 
oscillations. The example of this oscillations type are 
i^fj.,T or Vt conversions that are the possi- 

ble explanations of the solar and atmospheric neutrino 
problems. The second type of neutrino oscillations is the 
transitions between neutrino helicity states within one 
fiavor, or neutrino spin oscillations. The last opportunity 
is the combination of the two first neutrino oscillations 
types, i.e. the transition when both fiavor and helicity 
states can change. This type of neutrino oscillations is 
called neutrino spin-flavor oscillations. It can be impor- 
tant in solving of the solar neutrino problem. In this 
paper we restrict ourselves to studying of neutrino spin 
oscillations. 

It was realized more than twenty years ago that ex- 
ternal fields drastically change the process of vacuum 
neutrino oscillations. In Refs. 0, Q it was established 
that neutrino weak interactions with background mat- 
ter result in the resonant enhancement of neutrino fia- 
vor oscillations. The electromagnetic interactions are of 
great importance in studying of neutrino spin and spin- 
fiavor oscillations. For example, the neutrino interaction 
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with an external electromagnetic field provides one of 
the mechanisms for the mixing between different helicity 
eigenstates. It is also interesting to study the infiuence 
of gravitational fields on neutrino oscillations (neutrino 
does not seem to participate in strong interactions). De- 
spite the gravitational interaction is relatively weak com- 
pared to electromagnetic and weak interactions there are 
rather strong gravitational fields in the Universe. Thus 
neutrino oscillations in gravitational fields are of interest 
in astroparticle physics and cosmology. 

The infiuence of the gravitational interaction on neu- 
trino oscillations has been studied in many publications 
(see, e.g., Refs. |llilllili|l3,|llli3)- First we men- 
tion Refs. 0, where the formalism for the description 
of neutrino oscillations in gravitational fields in presence 
of matter and external electromagnetic fields was worked 
out. One of the approaches to the description of the spin 
dynamics of a neutrino in a gravitational field consists 
in decomposing of the Dirac Hamiltonian (in presence 
of an external gravitational field) and establishing terms 
that mix different chirality components of the neutrino 
wave function. This method was used in Ref. To 
describe the helicity evolution in a gravitational field one 
can use the evolution equation in Heisenberg representa- 
tion with the Hamiltonian in Foldy-Wouthuysen repre- 
sentation. This approach was implemented in Ref. 
Neutrino spin oscillations and spin light of a neutrino in 
weak gravitational fields were studied in our work |l2j |. 
Gravity induced neutrino fiavor oscillations were exam- 
ined in Refs. HH. 

There is, however, another approach for the treatment 
of the spinning particle dynamics in an external gravita- 
tional field. The system of equations for the description 
of the spinning particle motion in gravitational fields was 
derived by A. Papapetrou in Ref. 13] . Since then plenty 
of works where Papapetrou equations were analysed have 
been published. The major difficulties in solving Papa- 
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petrou equations are that these equations are non-linear 
in the particle spin as well as the motion of a spinning 
particle deviates from the geodesical. Thus to describe 
the motion of a finite size spinning particle in a gravita- 
tional field one has to take into account the contribution 
of the spin term into the particle motion law. It was es- 
tablished in Ref. that Papapetrou equations can be 
significantly simplified when one uses the quasi-classical 
approach and linear in spin approximation of these equa- 
tions for a point particle. 

In this paper we study neutrino spin oscillations in 
gravitational fields within the quasi-classical approach. 
Note that the quasi-classical treatment of neutrino spin 
oscillations in moving and polarized matter under the 
influence of electromagnetic fields was considered in 
Refs. Hi IB 113. In Sec. ^ we formulate the main equa- 
tions necessary for the description of the neutrino motion 
and spin evolution. The limit of a weak gravitational 
field is examined in Sec. IIIII On the basis of the gen- 
eral equation presented in Sec.|n]we obtain the effective 
Hamiltonian for neutrino spin oscillations. This effective 
Hamiltonian is valid for arbitrary neutrino velocity. We 
also compare the neutrino spin evolution equation with 
the analogous one derived in our previous work [T^ and 
find out that the equation obtained in the present work 
correctly accounts for the neutrino velocity dependence. 
Then we receive neutrino transition probability when a 
particle propagates in the gravitational field of a rotating 
massive object. In Sec. IIVI we apply the general tech- 
nique to the description of neutrino spin oscillations in 
the Schwarzschild metric. We obtain the neutrino spin 
evolution equation which is valid for the neutrino motion 
even in the vicinity of a black hole. The effective Hamilto- 
nian and the transition probability are also derived. We 
examine neutrino oscillations process on different circular 
orbits and analyze the frequencies of the spin transitions. 
In Sec. El we discuss our results. The calculation of co- 
variant derivatives of the vierbein vectors f Appendix 
and the basic elements of SL{2,C) group f Appendix IB|) 
are also presented. 



II. MOTION OF A SPINNING POINT 
PARTICLE IN GRAVITATIONAL FIELDS 

The motion of a spinning particle in gravitational fields 
was described in Ref. 13] . The evolution of the spin 
tensor of a particle 5'"^ and particle's momentum is 
described by the following equations, 
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where u'' is the unit tangent vector to the center-of-mass 
world line, r is the parameter (not necessarily the proper 
time) which changes along the world line, D/Dt denotes 
the covariant derivative along the world line and R'^^ pa- 



is the Riemann tensor. The spin vector is introduced in 
the following way (see, e.g., Ref. 0|), 
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where s^i,xp is the completely antisymmetric tensor den- 
sity, g = det(5^^) and — p^p'^ is the constant of the 
particle's motion. 

When we study the motion of point particles in a grav- 
itational field, Eqs. H2.1|l and ()2.2ll can be significantly 
simplified. It was shown in Ref. [ISj that using only the 
Principle of General Covariance the equations for the de- 
scription of particle's spin vector S'^ and four-velocity 
= dx^/dr evolution take the form. 
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Here the covariant derivatives are taken with respect to 
the proper time r. It follows from Eqs. H2.3|l and H2.4|l 
that particle's spin and four-velocity are parallel trans- 
ported along its world-line. It should be noted that usual 
properties of the spin vector in Minkowski space-time, 
namely S^Ufj_ — and S'^S'^ is the constant value, remain 
unchanged when a particle moves in curved space-time. 
One can verify these properties using Eqs. H2.3|l and (|2.4|l . 
The basic Eqs. (gS and (E31l can be also rewritten with 
help of the Christoffel symbols of the second kind T^afi, 
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Eq. (|2.5|) describes the spin evolution in a general coor- 
dinate system. However the particle's properties are de- 
termined by the spin vector components measured with 
respect to the particle's rest frame. The particle's rest- 
frame-spin precession in gravitational fields was discussed 
in Ref. .141]. In order to proceed in our analysis of par- 
ticle's spin evolution, we should rewrite Eqs. (|2.5|l and 
(|2.t)|) in a locally minkowskian frame. One can imple- 
ment this coordinates transformation with help of the 
vierbein vectors e"^. We recall the basic vierbein vectors 
properties. 



(2.7) 
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where Cq^ = Vabg^'^e u is the inverse vierbein and rjab = 
diag(-|-l, — 1, — 1, — 1) is the metric tensor in a locally 
minkowskian frame. 

The components of the spin and four-velocity with re- 
spect to a locally minkowskian frame are. 
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The equations for the description of the s° and u" evo- 
lution have the form (see also Ref. (141). 



III. NEUTRINO SPIN EVOLUTION IN WEAK 
GRAVITATIONAL FIELDS 



where G"'' = r]'"'Tj'"^ 



lu ~ 
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~dt ~ 
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= — G''" is the analog 


; of the 



electromagnetic field tensor and ^abc 
are the Ricii rotation coefficients. In Eqs. H2.8|l and (|2.9|) 
we keep the usual notation ^ — — dt/dr, however it 
should be noted that 7 7^ The derivatives in the right- 
handed sides of these equations are taken with respect to 
the laboratory time. 

To make the coordinate transformation into the par- 
ticle's rest frame one uses a boost within a locally 
minkowskian frame. The relation between the spin vec- 
tor s° and the three dimensional spin in the particle's 
rest frame is given by following expression, 



(Cu),C 



u(Cu) 

1 + 
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Here it*^ and u are the time and space components of the 
four- velocity u". We remind that u° is the four velocity 
vector in the vierbein frame. 

Using Eqs. and ifTTUIl we can readily derive the 

equation for the description of the three dimensional spin 
vector evolution. 



where 
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In deriving of Eqs. (|2.11|) and (|2.12|) we use the 
fact that any antisymmetric tensor in four dimensional 
(minkowskian) space-time can be expressed in terms of 
the two three dimensional vectors (analogs of electric and 
magnetic fields), i.e. Gab ~ (E,B), where Goi — Ei and 
Gij = —EijkBk- Eqs. H2.11|l and (|2.12l) describe particle's 
spin precession in arbitrary gravitational field. These 
equations are linear in spin vector. One can, however, 
discuss further (non-linear in spin) terms contributing to 
the particle's spin evolution (see Ref. [3). 

At the end of this section me mention that Eqs. (|2.8|) 
and (|2.9|l [or alternatively Eqs. H2.11|l and (|2.12|) ] are sim- 
ilar to the equations for the description of the spin and 
four-velocity evolution of a charged particle with gyro- 
magnetic ratio g = 2 interacting with an external electro- 
magnetic field (see also Ref. |l^). The spin evolution of a 
charged particle with g = 2 was studied in Refs. [20ll2]| . 
It was shown that the spin dynamics was completely de- 
termined by the particle's motion law, i.e. by the solution 
of the Lorentz equation (|2.9|) . 



In this section we study the particle spin precession in 
a weak gravitational field. We describe the neutrino spin 
oscillations and discuss the obtained results as well as we 
consider analogous approach elaborated in our previous 
work. 

Weak gravitational fields can be found, for instance, at 
great distances from a finite massive object under study. 
In this case we can always choose the quasi-minkowskian 
coordinate system. 



(3.1) 



where 77^^ = diag(+l, —1, —1, —1) is the Minkowski met- 
ric tensor. The metric perturbation hf^, has to vanish at 
the infinity. One of the possible examples of the metric 
given in Eq. H3.1|) is the gravitational field created by a 
massive rotating object at great distances. Using post- 
Neutonian approximation we obtain for the components 
of the tensor h„,,y (see, e.g., Ref. l2^). 
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i.f^i, (see, e.g., Ref. |22 

= 2(/3, h,j = 2ipS,j, 



(3.2) 
(3.3) 



where M is the mass of the object, J is its total angular 
momentum and r is the distance from the object. 

First let us consider the particle's spin evolution when 
an observer is in the particle's rest frame. For this pur- 
pose we can use directly Eq. (|2.5|l . In our case the four 
dimensional spin vector is reduced to the three dimen- 
sional one, i.e. (0, C)- Eq. (|2.5|l . written for the 
spatial components of the spin vector, takes the following 
form. 



(3.4) 



Here we assume that 17° = 1 and U = in particle's rest 
frame. Note that we do not distinguish between upper 
and lower indexes when we use three dimensional vectors. 
Christoffel symbols F^oj can be calculated with help of 
Eq. 



dhj 



dhj_ 

dx. 



= -o£yfe[Vxh]fe. (3.5) 



Using Eqs. H3.4I) and (|3.5|) we obtain the equation for the 
particle's spin evolution in its rest frame. 



dc 

dt 



-[Cx [Vxh]]. 



(3.6) 



One can see that Eq. 1)3. 6|) coincides (to within the sign) 
with the analogous equation derived in our previous work 
(see Ref. 12]) if we set there 7=1 and (3 — {). 

Now let us discuss the particle's spin precession in a 
weak gravitational field, with an observer being placed 
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in the laboratory frame, i.e. a particle has the non-zero 
velocity with respect to him. In this case we should use 
Eqs. (|2.11|l and (|2.12(l to describe particle's spin evolu- 
tion. The Ricci rotation coefficients were calculated in 
Ref. il4i | in the weak gravitational field approximation. 
They can be expressed in the following way, 

1 / dhbc dhac\ ,„ 

where hab and a;" are the components of the tensor /i^i^ 
and vector a;^ in the vierbein frame. Note that we do not 
distinguish the vierbein and world indexes in the weak 
gravitational field approximation. 

Using Eqs. H3.2|l and H3.7|) as well as the definition of 
the tensor Gab we receive the expressions for the "elec- 
tric" and "magnetic" fields, 

E = -jV(p+^v,Vh„ B = 7[vx V]</7+|[Vxh]. 

One can readily find the vector G, which determines the 
particle's spin precession, in the explicit form, 

G=|Q[Vxh] + ^2±I[vx V]^ (3.8) 
I 



It can be seen that this expression for the vector G agrees 
with the analogous formula obtained in our prevoius work 
|12| only for slow neutrinos having v <^ 1 (see also above). 
This discrepancy results from the incorrect account of the 
spin-orbital interaction in Ref. |12| |. However, if we study 
only the radial neutrino motion, the main contribution 
to the neutrino spin evolution comes from the interaction 
with the angular momentum of the massive object and 
we reach an agreement with our recent work. 

We consider therefore a neutrino propagating along the 
radius, i.e. v = v{r/r). In this case we can rewrite 
Eq. H3.10|l in the form, 

G=^{ij-n(Jn)(2+i)}, (3.11) 

where n is the unit vector along the neutrino veloc- 
ity. One can see that the vector G given in Eq. (|3.11|) 
coincides with that derived in our previous work (see 
Ref. [13 )■ Thus all results concerning the neutrino spin 
light obtained in Ref. 0| are valid only for the radial 
neutrino motion. 



where v = r is the particle's velocity and 7 = (1 — 
w2)-i/2. Note that if we set v = in Eq. ipP^ . we 
obtain the expression consistent with Eq. (|3.6() obtained 
directly from Eq. (|2.5|) . It proves the validity of the used 
technique. 

Eqs. H2.11(l and 1)3. 8|1 govern particle's spin evolution in 
a weak gravitational field. These equations are valid for 
arbitrary particle velocities. Using Eqs. (|2.11() and H3.8() 
we can describe neutrino spin oscillations in the weak 
gravitational field. For example, it is possible to study 
the particle's spin dynamics of a neutrino emitted in the 
vicinity of a rotating black hole and then propagating 
faraway from the massive object (see Ref. U^). Suppos- 
ing that a neutrino has the velocity directed along the 
base vector we obtain from Eq. (|3.8ll the expression for 
the effective Hamiltonian (see Ref. |15|) 

JYcff = --(tG), (3.9) 



where cr — (ci, (T2, fa) are the Pauli matrices. 

Let us consider neutrino spin oscillations when a neu- 
trino interacts with the gravitational field having the 
properties given in Eq. (|3.3|l . Using Eqs. (|3.3II and 1)3. 8|l 
we receive for the vector G the following expression. 



+ ("0) 

7 + 1 r'^ J 
I 

Now let us write down the effective Hamiltonian for 
neutrino spin oscillations in our case. We can always 
choose the coordinate system so that, 

n= (0,0,1), J = (Ji,0,J3), (3.12) 
Ji = Jsin-d, J3 ~ Jcos-d, 

where is the angle between vectors n and J. Using 
Eqs. H3.9|l . H3.11|l and 1)3.12(1 we obtain the expression for 
the neutrino spin oscillations effective Hamiltonian, 

i/e. = 4 f . r,f, , ''^''/^'^A . (3.13) 

7-3 \^smw/(27) — cosw J ^ ' 

Here r = Tq + vt, where Tq is the initial neutrino coordi- 
nate. It follows from Eq. 1)3. 13(1 that transitions between 
different polarization states of a neutrino are suppressed 
by the factor 1 /7 which is small for ultrarelativistic neu- 
trinos. The most intensive neutrino spin oscillations take 
place when a neutrino is propagating in the direction per- 
pendicular to the vector J (i? = 7r/2). There are no os- 
cillations at all when a neutrino is propagating along the 
vector J (i? = 0). 
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l[z;2j-v(vJ)] - A[vxr](v[rx J]) 
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The transition probability for the vl ^ oscillations 
can be obtained with help of the effective Hamiltonian 
given in Eq. Let us suppose that vr{Q) = and 

'^l(O) = 1, then we receive the expression for the prob- 
ability to find a right-handed neutrino in this neutrino 
beam, 



P{t) 



SVC? 1? 



cos^ d + sin i9 



sin^[$(r) 



where 



$(r) 



J /I 1 



(3.14) 



(3.15) 



X y 472 cos^ 'd + sin d. 



Here rg is the distance between a neutrino and the center 
of the massive object at t = 0. Eqs. (|3.14|) and H3.15|l co- 
incide with the similar expression for the transition prob- 
ability derived in our previous work [12j |. 



IV. NEUTRINO SPIN OSCILLATIONS IN 
SCHWARZSCHILD METRIC 



satisfy the general properties of the vierbein vectors given 
in Eq. 

Using Eqs. we can find vierbein components 

of the vector u", 

u'^ = {-iA,UrA-^,Uer,U4,rs\xv9). 

With help of Eqs. as weU as using the defini- 

tion of the tensor Gab one can write down the components 
of the "electric" and "magnetic" fields, 



E = 

B = {U^ cos 



(4.8) 



'U^Asine,U0A). 



Neutrino spin precession is determined by the vector n = 
G/7. The components of this vector can be found on the 
basis of Eqs. and ifl^H)) . 



r^i = — W0 cos 0, 



(4.9) 
(4.10) 

(4.11) 



General equations derived in Sec. allow one to de- 
scribe particle's spin evolution not only in weak gravita- 
tional fields. In this section we discuss the neutrino spin 
evolution and oscillations when a particle interacts with 
a non-rotating black hole in the case when a neutrino 
propagates even in the vicinity of a black hole. 

When we study the gravitational field of a non-rotating 
black hole, the interval is known to be expressed with 
help of the Schwarzschild metric. 



where 



A^dt^ 



A-^dr^ 



A^ 



\de^ 



(4.1) 



and rg is the Schwarzschild radius. In Eq. (|4.1(l we use 
the spherical coordinate system for the world (not for 
the vierbein) coordinates. We can decompose the four- 
velocity [/^ with help of the spherical coordinates. 



- {U\Ur,Ug,U^). 



(4.3) 



In order to construct the "electric" and "magnetic" 
fields we should choose the appropriate vierbein vectors. 
One can verify that the following vectors, 



(A 0,0,0), 

(0,^"\0,0), 
(0,0, r,0), 
(0,0,0,rsin6l), 



(4.4) 
(4.5) 
(4.6) 
(4.7) 



where v = (vr,vg,V(j,) are the components of the world 
(not the vierbein) velocity. In deriving of Eqs. H4.9(l - 
(14.11(1 we use the fact that U = 7V. 

Let us discuss a neutrino orbiting a black hole. For 
simplicity we consider only circular orbits with the ra- 
dius R. We may restrict ourselves to the studying of 
the orbits lying only in the equatorial plane (0 = 7r/2 or 
equivalently vg = 0) because of the spherical symmetry 
of the gravitational field. In this case fii = and = 
in Eqs. H4.9|l and H4.11|l . The expressions for the neutrino 
angular velocity and 7 are presented in Ref. p3 |. 



(4.2) and 



7 



dr 



1 



3rg 
2R' 



(4.12) 



(4.13) 



It should be noted that the vierbein four velocity now 
takes the form, m° — (7^, 0, 0, 7f0r). One can ver- 
ify that = 1 with help of Eqs. 1(4.21) . 1(4.12(1 and 
((4.13() . We also mention that Eq. ((2.9(1 is also identi- 
cally satisfied. Indeed GabU^ = since EivP = B2U^ and 
hence du°-/dT = 0. Therefore a neutrino has constant 
four velocity with respect to the vierbein frame. Using 
Eqs. ((4.21) . ((4.1011 and 114.13(1 we can rewrite in the 
more simple form. 



' 2 ■ 



(4.14) 



This equation has very clear physical meaning. It is well 
known that the spin precession of a charged particle with 
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g = 2 in the external electromagnetic field is completely 
determined by the particle's motion law (see Sec. ^ or 
Ref. 20]). Therefore Eq. H4.14|l is nothing else as the 
Lorentz transformation of the spin rotation frequency 
from the rest frame to the laboratory frame. 

We suppose that initially a neutrino is left-handed, i.e. 
its initial spin vector is antiparallel to the particle's ve- 
locity. According to Eqs. (|2.11|) and H4.9|l - H4.11|l the neu- 
trino spin rotates around the second axis. Therefore us- 
ing Eq. (|2.11|l we can construct the effective Hamiltonian 
for the neutrino spin oscillations in the gravitational field 
of a non-rotating black hole, 



0.07 



cff 



-iVL2 
iVL2 



(4.15) 



It is interesting to note that in Eq. H4.15|l the parameter 
|f72|7^0 if 7^ ^7^0 (see below). One can see it directly in 
Eq. 14.14|l . Using Eq. H4.15|l we write the expression for 
the neutrino transition probability. 



P{t) = sin^{n2t). 



(4.16) 



We can see that there is a full mixing in our case and 
neutrino transition probability can achieve a unit value. 

It is possible to describe the neutrino spin evolution 
using formalism elaborated in Ref. ^21,] where the the 
approach for the studying of Lorentz and Bargmann- 
Michel-Telegdi equations was described. It follows from 
the results of that work that the spin dynamics is com- 
pletely determined by the particle's motion law when 
g = 2. Eq. H2.9|l can be rewritten in the equivalent form, 



dT 



A = FA, -cr(B - iE) 



and the operator A implements the shift of the particle's 
velocity along the trajectory. We present the main prop- 
erties of the SL{2, C) group in Appendix IbI If one has 
found the operator A, the resolvent of the Eq. H2.8|l takes 
the form (g = 2), 



R{t,to)^ L-'{t)AL{to), 



where the matrix L implements the boost from the rest 
frame to the laboratory frame. 

If we consider the neutrino motion in the equatorial 
plane from the very beginning, the operator A is ex- 
pressed in the following way. 



A = exp(FT) = ( ''"''."^ '^"""^^ 
~ ' la_sinat cosai j 



(4.18) 



where a 



VW~W/{2-f) and 

E±B 



Note that using Eqs. and (jlT^ we can 

show that a — Vtp/{2^). It can be verified that AuA^ — u, 
where 

' + 




FIG. 1: Neutrino spin oscillations frequency versus the radius 
of the neutrino orbit. 



i.e. a neutrino moves along a strait line (in the veirbein 
frame). This fact also agrees with the result obtained 
above in present paper. The operator L does not depend 
on time since the neutrino's four velocity is constant. 
This matrix is expressed in the following way (see Ap- 
pendix 0, 



L = 







(4.19) 



Using Eqs. (|4.17() - (|4.19() we receive the resolvent R{t) in 
the form (here we assume that to = 0), 



Bit) 



cos at sin at 
— sin at cos at 



(4.20) 



In order to find the spin rotation axis one adopts the unit 
vector k along this axis, so that [R, crk] = 0. Here use 
the fact that R G SO (3) and = R~^ ■ It can be easily 



(4.17) derived with help of Eq. (|On|) that [R,ay\ = 0. Thus 



the spin precession axis is the second axis in the complete 
agreement with the above obtained results. 

One can also plot the frequency of neutrino spin os- 
cillations versus the radius of the orbit. It is possible to 



see on Fig. [T]that |f72| = at R = l.brg 



at R 



oo. One can conclude that IfJo 



and \VL2\^ Q 
has its maxi- 
= 2rg [see also 



mal value (equal to 6.25 x 10~^r~^) at R 
Eq. 14.1411 ]. Let us evaluate the number of revolutions 
N , that a neutrino should make, necessary for the total 
spin flip. Using Eqs. H4.12|l and (|4.14|l we can easily find 
for iV, 



N = 



V4, 



i\r- 



= 2. 



It is interesting to evaluate the characteristic period of 
the neutrino spin oscillations (oscillations length). For 
M = IOMq at i? = 2rg we get for T = Tr/lf^a] ~ 4.94 x 
10-3 s. 
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We can discuss another orbit, namely R = l.Sr^. It fol- 
lows from Eq. H4.14|l (see also Fig. that the frequency 
of the transitions between two polarization states van- 
ishes. It should be noted that dr — a,t R — l.5rg [see 
Eq. H4.13|l ]. Therefore this orbit corresponds to mass- 
less particles. We reveal that the property of a massless 
particle in the Minkowsky space-time to keep unchanged 
its spin direction with respect to the momentum remains 
unaffected in the Schwarzschild metric. It is interesting 
to note that the absence of the spin rotation for massless 
particles is valid only for the Schwarzschild metric. It 
was shown in Refs. |24Ll25| that for a background metric 
with at least one off-diagonal space component, e.g. the 
Kerr metric, in the limit — )■ still there are tran- 
sitions between neutrino and antineutrino, i.e. neutrino 
spin flip. This phenomenon is used in Refs. ,24> .2&] to 
explain the neutrino- antineutrino asymmetry problem. 

It is also worth examining the behaviour of for large 
orbits. If a neutrino moves on a remote orbit, the gravi- 
tational field is weak. Thus we again do not distinguish 
between world and vierbein indexes. We define the co- 
ordinate system so that the angular velocity is expressed 
in the following way, Wfi-amo = v^ez , where is the unit 
vector. Using Eq. (|4.14|l in the weak field limit we obtain 
for the spin rotation frequency, 



it for our purposes) is satisfied, 



'^abs 



1 



3rg 
4 R 



since the second vierbein coordinate axis is nothing else 
as the 0-axis and for the orbit in the equatorial plane we 
have eg = — e^. Note that the effective Hamiltonian in 
Eq. 1)4. 15|) determines the evolution of the neutrino rest- 
frame-spin vector. However an observer is in the labo- 
ratory frame {t is time in the laboratory frame). Thus 
f22 is the oscillations frequency measured in the labo- 
ratory frame. In order to obtain the spin rotation fre- 
quency measured by an observer in the co-moving frame 
we should subtract Wframc from Wabs since Wframc is the 
angular velocity of the orbital motion, i.e. the angular 
velocity of the rest frame rotation with respect to the 
laboratory frame. Finally we get for = i^ahs — '^^tramc 
the following expression, 



3rg 
4 R 



which agrees with the classical result for the spin rotation 
frequency in the co- moving frame (see, e.g., Ref. [glj or 
Eq. ifTTUIl in the hmit 7-^1). 

At the end of this section we briefly examine the ap- 
plicability of the quasi-classical approach to the descrip- 
tion of the neutrino spin evolution in the gravitational 
field of a black hole. Basing on the similarity between 
particle spin precession in gravitational and electromag- 
netic fields we can apply the results of Ref. where 
the quasi-classical approximation for the electron spin 
precession in an external magnetic field was considered. 
It was found in that paper that one can use the quasi- 
classical limit if the corresponding inequality (we rewrite 



h 

2£ 



dC 

dt' 



< 1 



(4.21) 



where £ is the neutrino energy in the vierbein frame, the 
derivative is also taken with respect to the vierbein time 
t'. We again study the case of circular orbits. For the 
neutrino energy one finds, £ = rrii^vP = m^7A, where 
mu is the neutrino mass. Using Eq. (|4.4|l we obtain that 
dt' — Adt. Let us discuss an electron neutrino with mass 
r7iy_. « 2 eV rotating around a black hole with M — \OMq 
on the orbit with radius R ~ 2rg. Then Eq. H4.21|l reads, 



h_ 

2£ 



dC 
dt' 



We can see that the condition of the validity of the quasi- 
classical approach is satisfied almost in all reasonable as- 
trophysical objects. 



V. CONCLUSION 

In conclusion we note that neutrino spin oscillations 
in gravitational fields within the quasi-classical approach 
have been studied. We have started from the particle 
spin evolution equation in the vierbein frame. On the 
the basis of this approach we have investigated neutrino 
spin oscillations in a weak gravitational field created by a 
massive rotating object. The three dimensional spin evo- 
lution equation has been obtained. We have improved 
analogous calculation performed in our previous work 
|12|. The neutrino spin evolution equation derived in 
the present work is valid for arbitrary neutrino velocities 
and correctly takes into account the neutrino velocity 
dependence. Then we have written down the neutrino 
oscillations effective Hamiltonian [Eq. I|3.13(l ] and transi- 
tion probability [Eqs. (|3.14(l and (|3.15() ] when a neutrino 
propagates along the radial direction. The neutrino spin 
conversion rate has been analyzed for different neutrino 
velocity directions. Then we have examined neutrino spin 
oscillations in the gravitational field of a non-rotating 
black hole. We have studied a neutrino moving on circu- 
lar orbits, with our calculations being valid for the neu- 
trino motion in the vicinity of a black hole. We have 
derived the effective Hamiltonian [Eq. H4.15|) ] and tran- 
sition probability [Eq. H4.16|l ] in this case. It has been 
demonstrated that neutrino spin oscillations occur in the 
Schwarzschield metric. We have also analyzed the de- 
pendence of the neutrino spin oscillations frequency on 
the radius of the orbit. The validity of the quasi-classical 
approach to the description of neutrino spin oscillations 
in a gravitational field has been analyzed in our work. 

It should be noted that the method used in the present 
paper for the studying of neutrino spin oscillations is 
valid not only for neutrinos but also for any spin- 1/2 
particles. The basic equations (|2.11|) and H2.12|l were 
demonstrated in Ref. |22] to be applicable to a massive 
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Dirac particle in a weak static gravitational field. There- 
fore we can study, e.g., electron spin evolution with help 
of this approach. It is known that rather frequently mas- 
sive gravitational objects have inherent magnetic fields 
and are surrounded by the moving matter (accretion 
disks). Thus, if we discuss the particle spin oscillations 
process comprehensively, i.e. take into account all fac- 
tors, the cases of neutrinos and electrons are different 
because these particles have diverse types of the interac- 
tion with magnetic field and with the background matter. 
In principle massless particles can be also included into 
consideration (see Sec. IIV|) . Despite of the fact that the 
mass of a particle is absent in Eqs. I|2.11|l and H2.12|l . 
the law of motion (i.e. the function u°-{r,t), which is 
involved in the spin evolution equations) of massless par- 
ticles differs from that of massive particles. However in 
the case of massless particles the helicity and chirality 
are the same and one can apply the methods used in the 
recent work |26l |. 

According to Eqs. H4.9|l - I|4.11|l there is no neutrino 
spin flip in case of radial neutrino motion in the 
Schwarzschield metric. This our result agrees with the 
deduction of Ref. "T] where the radial propagation of neu- 
trinos from active galactic nuclei was studied. However 
one can hardly agree with the claims of Ref. 01 that 
a spherically symmetric, static Schwarzschild space-time 
cannot cause the particle spin flip. That statement con- 
front both our results and the analysis of Ref. 29] in 
which it was shown by means of direct calculations that 
in the Schwarzschield metric the spin of a Dirac fermion 
can precess and change its direction with respect to the 
particle momentum. The calculations in that paper were 
based on Dirac equation in curved space. That result also 
agrees with the statement of our paper that the gravita- 
tional field of a non-rotating black hole can cause neu- 
trino spin oscillations. 

Note that the results of our paper can be applied to the 



description of the neutrino spin evolution in various as- 
trophysical media. It is known that relic massive neutri- 
nos can cluster into a halo around a galaxy contributing 
to cold dark matter. Despite of the fact that the average 
size of a halo ranges between 10^ and 10'^ k pc for a typical 
galaxy like the Milky Way (see, e.g., Ref. |33|), there is a 
possibility for a neutrino to be gravitationally captured 
on orbits close to the central massive object. In this case 
strong gravitational field can essentially contribute to the 
neutrino spin oscillations process. 
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APPENDIX A: COVARIANT DERIVATIVES OF 
THE VIERBEIN VECTORS IN 
SCHWARZSCHILD METRIC 

The covariant derivatives of the veirbein vectors enter 
in the definition of the tensor Gat- In this appendix we 
present the main formulae necessary for the covariant 
derivatives calculation. The covariant derivative of e^^ 
is equal to 

^afi;u — qJ^ — r ^i/BaA, (Al) 

where ^i, are the Christoffel symbols. The non-zero 
Christoffel symbols for the Schwarzschield metric are 
(see, e.g., Ref. HJ), 



yr _ • g 
r-^^e = cot0. 



1 



r) 

— 

-L fAfA — 



' sm y cos ( 



2 n T^t 



r''00 = [vg - r) sin 6, T tr 



2r{r 



(A2) 



Using Eqs. (|4.4() - (|4.7|) . ((XlT) and (|X2|) we obtain the ex- 
pressions for the covariant derivatives in the following 
form, 



0\ 

-{rg/2r^)A-^ 



0/ 



(A3) 



/(rg/2r2)A 





{rg - r)^-i 

0- 

^0 

1 



VO —rsinOcosBj 



(A4) 



(A5) 
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/O 

sin^ 

rcos0 

\0 



(A6) 



where the first index (/i) numbers lines and the second 
one (v) numbers the rows of the matrix. 

It is convenient to rewrite the definition of the tensor 
Gab in the following way, 



Gab 



;ueb^Ur 



(A7) 



Here we used the properties of the vierbein vectors [see 
Eq. 1)2.71) ]. Eq. (|A7|) is more convenient for the further 
calculations since it allows one to express the tensor Gab 
directly in terms of the components of the four vector 
rather than 



(det A = 1 and A°o > 1 for A € In this appendix 

we briefly discuss the local isomorphism between these 
groups (see also Ref. js^)- 



For any element x'^ of the minkowskian space we set 
the matrix x = x'^e^, where eo is the 2x2 unit matrix 
and e = (T. The Lorentz transformation acts on x as 
Ax — AxA^, where A^ is the Hermitian cojugate to the 
matrix A. For instance, the Lorentz transformation of 
the form, 



ip =x'^ coshx + (xl) sinhx, 
y =x — (xl)l + [(xl) coshx + x'^ sinhxjl, 



APPENDIX B: UNIVERSAL COVERING OF THE 
LORENTZ GROUP 

It is well known that SL{2, C) is the universal covering 
group for the orthochronous proper Lorentz subgroup L \ 



corresponds to the matrix H{\,x) = exp[x(crl)/2]. Here 
1 is the three dimensional unit vector and x is the real 
parameter. 
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